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Aims. We probe the dependence of a 2 /fi on the ambient matter density by means of submm- and mm-wave bands 
^"P" spectral observations in the Milky Way. 

Methods. A procedure is suggested to explore the value of F = a 2 /fi, where /i = m c /m p is the electron-to-proton 
. — , mass ratio, and a = e 2 /(he) is the fine-structure constant. The fundamental physical constants, which are measured 

in different physical environments of high (terrestrial) and low (interstellar) densities of baryonic matter are supposed 
to vary in chameleon-like scalar field models, which predict that both masses and coupling constant may depend 

Oon the local matter density. The parameter AF/F — (Jobs ~ -Fiab)/-Fiab can be estimated from the radial velocity 
offset, AV = Kot — Vf s , between the low-laying rotational transitions in carbon monoxide 13 CO and the fine-structure 
\,J ■ transitions in atomic carbon [Cl]. A model-dependent constraint on Aa/a can be obtained from AF/F using Afi/fi 

independently measured from the ammonia method. 

Results. Currently available radio astronomical datasets provide an upper limit on \AV\ < 110 m s _1 (l<r). When 
interpreted in terms of the spatial variation of F, this gives \ AF/F\ < 3.7 x 10~ 7 . An order of magnitude improvement 
of this limit will allow us to test independently a non-zero value of Afi/fi — (2.2±0.4 sta t ±0.3 sys ) x 10 -8 recently found 
with the ammonia method. Taking into account that the ammonia method restricts the spatial variation of p at the 
level of \A/i/n\ < 3 x 10~ 8 and assuming that AF/F is the same in the entire interstellar medium, one obtains that the 
spatial variation of a does not exceed the value |A«/a| < 2 x 10~ 7 . Since extragalactic gas clouds have densities similar 
to those in the interstellar medium, the bound on Aa/a is also expected to be less than 2 x 10~ 7 at high redshift if no 
significant temporal dependence of a is present. 

■ Key words. Line: profiles - ISM: molecules - Radio lines: ISM - Techniques: radial velocities 
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| 1. Introduction fission reactor which allows us to probe a(t) at At ~ 2 x 10 

. • yr (z ~ 0.4). Assuming possible changes only in the electro- 

^j- . The dimensionless physical constants like the electron-to- magnetic coupling constant, Gould et al. (2006) obtained 

O i proton mass ratio, fx = m c /m p , or the fine-structure con- a mo del dependent constraint on \Aa/a\ < 2 x 10~ 8 . 

stant, a = e 2 /(hc), are expected to be dynamical quanti- However, when the strength of the strong interaction, - 

ties in modern extensions of the standard model of parti- the parameter A QC d, ~ is also suggested to be variable, 

cle physics (Uzan 2003; Garcia-Berro et al. 2007; Martins the Oklo data does not provide any bound on the variation 

2008; Kanekar 2008; Chin et al. 2009). Exploring these pre- f a (Flambaum & Shuryak 2002; Chin et al. 2009). 

dictions is a subject of many high precision measurements Current astrophysical measurements at higher redshifts 

^ , in contemporary laboratory and astrophysical experiments. are as follows. There was a claim for a variability of a at the 

. The most accurate laboratory constraints on temporal a- 5<r confidence level: Aa/a = -5.7±1.1 ppm (Murphy et al. 

and ^-variations of a/a = (-1.6 ± 2.3) x 10~ 17 yr" 1 , and 20040, but this was not confirmed in other measurements 

(iff! = (1.6±1.7) x 10 15 yr 1 were obtained by Rosenband wn i c h led to the upper bound \Aa/a\ < 2 ppm (Quast et 

et al. (2008), and Blatt et al. (2008), respectively. a l. 2004; Levshakov et al. 2005; Srianand et al. 2008; Molaro 

In case of monotonic dependence of a(t) and /z(i) on cos- et al. 2008a). 

mic time, at redshift z ~ 2 (corresponding look-back time Measurements of the cosmological ^-variation exhibit a 

is At ~ 10 10 yr) the changes of a and n would be restricted similar tendency. Non-zero values of Ay/y = -30.5 ± 7.5 

at the level of \Aa/a\ < 4 x 10" 7 and |A/V/x| < 3 x 10~ 5 . ppm, Au./u. = -16.5 ± 7.4 ppm (Ivanchik et al. 2005), 

Here Aa/a (or A/if/i) is a fractional change in a be- and A\x/\x = -24 ± 6 ppm (Reinhold et al. 2006) found 

tween a reference value «i and a given measurement a 2 at z = 2.595 (Q 0405-443) and z = 3.025 (Q 0347- 

obtained at different epochs or at different spatial coordi- 383) from the Werner and Lyman bands of H 2 were later 

nates: Aa/a = (a 2 -a 1 )/a 1 . refuted by Wendt & Reimers (2008), King et al. (2008) 

These constraints are in line with geological measure- 

ments of relative isotopic abundances in the Oklo natural 1 Hereafter, 1 ppm = 10 -6 . 
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and Thompson et al. (2009) who used the same optical 
absorption-line spectra of quasars and restricted changes 
in /i at the level of A/1//1 < 6 ppm. The third H2 sys- 
tem at z = 2.059 towards the quasar J2123-0050 also 
does not show any evidence for cosmological variation in 
fi: Ajti/jU = — 5.6±5.5 sta t ± 2.9 sys ppm (Malec et al. 2010). 
More stringent constraints were obtained at lower redshifts 
from radio observations of the absorption lines of NH3 and 
other molecules: |A/x//x| < 1.8 ppm at z = 0.68 (Murphy et 
al. 2008), and |A/i//i| < 0.6 ppm at z = 0.89 (Henkel et al. 
2009). Two cool gas absorbers at z = 1.36 (Q 2337-011) and 
z = 1.56 (Q 0458-020) were recently studied in the H 1 21cm 
and C IAA1560, 1657 absorption lines providing a constraint 
on the variation of the product X = g p a 2 n (here g p is the 
proton gyromagnetic ratio): AX/X — — 6.8±1.0 s tat±6.7 sys 
ppm (Kanekar et al. 2010). Thus, the most accurate as- 
tronomical estimates restrict cosmological variations of the 
fundamental physical constants at the level of ~l-2 ppm. 

The estimate of fractional changes in Aa/a and Afx/fj, 
by spectral methods is always a measurement of the rel- 
ative Doppler shifts between the line centers of different 
atoms/molecules and their comparison with corresponding 
laboratory values (SavcdofF 1956; Bahcall et al.1967; Wolfe 
et al. 1976; Dzuba 1999, 2002; Levshakov 2004; Kanekar & 
Chengalur 2004). To distinguish the line shifts due to radial 
motion of the object from those caused by the variability of 
constants, lines with different sensitivity coefficients, Q, to 
the variations of /i and/or a are to be usecQ- It is clear that 
the larger the difference |AQ| between two transitions, the 
higher the accuracy of such estimates. 

Optical and UV transitions in atoms, ions and molec- 
ular hydrogen H2 have similar sensitivity coefficients with 
|A<5| not exceeding 0.05 (Varshalovich & Levshakov 1993; 
Dzuba 1999, 2002; Porsev et al. 2007). For atomic spec- 
tra, the estimate of Aa/a is given in linear approximation 
(|Aa/a| < 1) by (e.g., Levshakov et al. 2006): 

Aa ^ (V 2 - Vi) _ AV 

a ~ 2c(Qi - Ca) ~ 2cAQ ' [ ' 

where Vi,V2 are the radial velocities of two atomic lines, 
and c is the speed of light. It was shown in Molaro et al. 
(2008b) that the limiting accuracy of the wavelength scale 
calibration for the VLT/UVES quasar spectra at any point 
within the whole optical domain is about 30 m s _1 , which 
corresponds to the limiting relative accuracy between two 
lines measured in different parts of the same spectrum of 
about 50 m s _1 . Taking into account that |AQ| ~ 0.05, it 
follows from Eq.fTJ) that the limiting accuracy of Aa/a is 
2 ppm, which is the utmost value that can be achieved in 
observations of extragalactic objects with present optical 
facilities. 

A considerably higher sensitivity to the variation of 
physical constants is observed in radio range. For exam- 
ple, van Veldhoven et al. (2004) first showed that the in- 
version frequency of the (J,K) = (1,1) level of the am- 
monia isotopologue 15 ND3 has the sensitivity coefficient 
Q M = 5.6. Compared to optical and UV transitions, the 
ammonia method proposed by Flambaum & Kozlov (2007) 
provides 35 times more sensitive estimate of A\x/\x from 
measurements of the radial velocity offset between the NH3 

2 Q is a dimensionless coefficient showing a relative change of 
the atomic transition frequency uii in response to a change of 
the physical constant F: Awi/ui = QiAF/F. 



(J,K) = (1,1) inversion transition at 23.7 GHz and low- 
lying rotational transitions of other molecules co-spatially 
distributed with NH3: 

^ « 0.289^ . (2) 

fi c 

The ammonia method was recently used to explore pos- 
sible spatial variation^ of physical constants from obser- 
vations of prestellar molecular cores in the Taurus giant 
molecular cloud (Levshakov et al. 2010, hereafter L10), the 
Perseus cloud, the Pipe Nebula, and Infrared dark clouds 
(Levshakov et al. 2008b; Molaro et al. 2009). In contrast 
to the mentioned above laboratory constraints on temporal 
variations, this method reveals a tentative spatial variation 
of Afi/n at the level of Afi/ fi = (2.2±0.4 stat ±0.3 sys ) x 10~ 8 
(L10). The corresponding conservative upper limit in this 
case is equal to \A/j,//j,\ < 3 x 10~ 8 . 

In the present paper we consider fractional changes of a 
combination of two constants a 2 and fi, F — a 2 / fi, which 
are estimated from the comparison of transition frequencies 
measured in different physical environments of high (ter- 
restrial) and low (interstellar) densities of baryonic matter. 
The idea behind this experiment is that some class of scalar 
field models — so-called chameleon-like fields — predict the 
dependence of both masses and coupling constant on the 
local matter density (Olive & Pospelov 2008; Upadhye et 
al. 2010). Chameleon- like scalar fields were introduced by 
Khoury & Weltman (2004a,b) and by Brax et al. (2004) to 
explain negative results on laboratory searches for the fifth 
force which should arise inevitably from couplings between 
scalar fields and standard model particles. The chameleon 
models assume that a light scalar field acquires both an 
effective potential and effective mass because of its cou- 
pling to matter that depends on the ambient matter den- 
sity. In this way, the chameleon scalar field may evade lo- 
cal tests of the equivalence principle and fifth force exper- 
iments since the range of the scalar-mediated fifth force 
for the terrestrial matter densities is too small to be de- 
tected. Similarly, laboratory tests with atomic clocks for co- 
variations are performed under conditions of constant local 
density and, hence, they are not sensitive to the presence 
of the chameleon scalar field (Upadhye et al. 2010). This is 
not the case for space-based tests, where the matter density 
is considerably lower, an effective mass of the scalar field 
is negligible, and an effective range for the scalar-mediated 
force is large. Light scalar fields are usually attributed to 
a negative pressure substance permeating the entire visible 
Universe and known as dark energy (Caldwell et al. 1998). 
This substance is thought to be responsible for a cosmic 
acceleration at low redshifts, z < 1 (Peebles & Rata 2003; 
Brax 2009). 

2. [Ci] and CO lines as probes of a 2 //j, 

The variations of the physical constants can be probed 
through atomic fine-structure (FS) and molecular rota- 
tional transitions (Levshakov et al. 2008a; Kozlov et al. 
2008). The corresponding lines are observed in submm- and 
mm- wavelength ranges. Along with a gain in sensitivity, the 
use of such transitions allows us to estimate constants at 

3 Hereafter, the term 'spatial variation' means a possible 
change in fi between its terrestrial and interstellar values. 



S. A. Levshakov et al.: Searching for spatial variations of a 2 / 'fi in the Milky Way 



3 



very high redshifts (z > 5) which are inaccessible to optical 
observations. 

Let us consider radial velocity offsets between molecu- 
lar rotational and atomic FS lines, AV = V T0 t — Vfs- The 
offset AV is related to the parameter AF/F as follows 
(Levshakov et al. 2008a): 

AF/F « 2Aa/a - Afi/fi = AV/c . (3) 

The velocity offset in Eq.® can be represented by the sum 
of two components 

AV = AVf + AV n , (4) 

where AVf is the shift due to F-variation, and AV n is the 
Doppler noise — a random component caused by possi- 
ble local offsets since transitions from different species may 
arise from slightly different parts of a gas cloud, at different 
radial velocities. 

The Doppler noise yields offsets which can either mimic 
or obliterate a real signal. Nevertheless, if these offsets arc 
of random nature, the signal AVf can be estimated statis- 
tically by averaging over a large data sample: 

{AV} = (AV f ), Var(AV) = Var(AV f ) + Var(AV n ) . (5) 

Here we assume that the noise component has zero mean 
and a finite variance. 

The Doppler noise component can be minimized if the 
chosen species are closely trace each other. An appropriate 
pair in our case is the atomic carbon FS transitions and ro- 
tational transitions of carbon monoxide 13 CO. The spatial 
distributions of 13 CO and [Ci] are known to be well corre- 
lated (Keene et al. 1985; Meixner & Tielens 1995; Spaans 
& van Dishoeck 1997; Ikeda et al. 2002; Papadopoulos et 
al. 2004). The carbon-bearing species C°, C + , and CO are 
observed in photodissociation regions (PDRs) - neutral re- 
gions where chemistry and heating are regulated by the 
far-UV photons (Hollenbach & Tielens 1999). The PDR is 
either the interface between the H II region and the molec- 
ular cloud or a neutral component of the diffuse interstellar 
medium (ISM). Far-UV photons (6.0 eV < hv < 13.6 eV) 
are produced by OB stars. Photons with energy greater 
than 11.1 eV dissociate CO into atomic carbon and oxy- 
gen. Since the C° ionization potential of 11.3 eV is quite 
close to the CO dissociation energy, neutral carbon can 
be quickly ionized. This suggests the chemical stratifica- 
tion of the PDR in the line C + /C°/CO with increasing 
depth from the surface of the PDR. Then, one can as- 
sume that in the outer envelopes of molecular clouds neu- 
tral carbon lies within a thin layer determined by the equi- 
librium between photoionization/recombination processes 
on the C + /C° side, and photodissociation/molecule forma- 
tion processes on the C°/CO side. However, observations 
(Keene et al. 1985; Zhang et al. 2001) do not support such 
a steady-state model which predicts that C° should arise 
only near the edges of molecular clouds. To explain the ob- 
served correlation between the spatial distributions of C° 
and CO, inhomogeneous PDRs with clumping molecular 
gas were suggested. The revealed ubiquity of the [C i] tran- 
sition 3 Pi — » 3 Po in molecular clouds is in agreement with 
clumpy PDR models (Meixner & Tielens 1995; Spaans et 
al. 1997; Papadopoulos et al. 2004). 

The ground state of the C° atom consists of the 3 -Pi,2,3 
triplet levels. The energies of the fine-structure excited lev- 
els relative to the ground state are £o,i = 24 K, and 



£■0,2 = 63 K, and the transition probabilities are Ai q = 
7.932 x 10~ 8 s" 1 , and A 2 ,i = 2.654 x 10~ 7 s" 1 (Silva & 
Viegas 2002). The excitation rates of the [Ci] J — 1 and 
,7 = 2 levels for collisions with H2 at T^ n ~ 20 K are 
30,1 ~ 9o,2 ~ 10~ 10 cmV 1 (Schroder et al. 1991). This 
implies that for the J — I and J = 2 levels the critical 
densities are 1000 cm~ 3 and 3000 cm~ 3 , respectively. The 
low- J rotational transitions of CO trace similar moderately 
dense [n ~ 10 3 cm~ 3 ) and cold (Tki n ~ 20 K) gas. It is not 
completely excluded, however, that some heterogeneity of 
spatial distributions of [Ci] and 13 CO may occur resulting 
in the radial velocity offsets. 

In the chameleon-like scalar field models for density- 
dependent /i(p) and Oi(p) the fractional changes in these 
constants arise from the shift in the expectation value of 
the scalar field between high and low density environments. 
Since the matter density in the interstellar clouds is ~ 10 16 
times lower than in terrestrial environments, whereas gas 
densities between the molecular clouds themselves are much 
smaller (n gas ~ 10 3 — 10 5 cm~ 3 ), all interstellar clouds can 
be considered as having similar physical conditions irrespec- 
tive of their location in space. This means that the noise 
component in Eq. ([5]) can be reduced by averaging over indi- 
vidual AF/F values obtained from an ensemble of clouds 
for which the measurements of both [Ci] and 13 CO lines 
are available. 

Equations @ and show that in order to estimate 
AF/F and A/i//i with a comparable relative error the un- 
certainty of the velocity offset in (J3]) must be ~3.5 times 
smaller than that in the ammonia method (~5 m s -1 , see 
L10). At the moment such data do not exist. Both labora- 
tory and astronomical measurements of the [C 1] frequencies 
have much larger uncertainties. For example, the rest fre- 
quencies of the [Ci] J = 1 - transition 492160.651(55) 
MHz (Yamamoto & Saito 1991) and J = 2 — 1 transition 
809341.97(5) MHz (Klein et al. 1998) are measured with the 
uncertainties of e v — 33.5 m s _1 and 18.5 m s _1 , respec- 
tively. For 13 CO the rest hyper-fine frequencies of low- J ro- 
tational transitions are known with good accuracy: i>i-q = 
110.201354280(37) GHz, and v 2 -i = 220.398684129(66) 
GHz, i.e., e v < 0.1 m s _1 (Cazzoli et al. 2004). Suggesting 
that the laboratory error e v ~ 34 m s^ 1 dominates over the 
errors from AV measurements, one obtains that the AF/F 
limiting accuracy is 0.1 ppm. To put in another words, if 
both species arise from the same volume elements and their 
radial velocities are known with a typical error of ~ 100 
m s _1 (e.g., Ikeda et al. 2002), then the mean AV can be 
estimated with a statistical error of ~ 30 m s _1 from an 
ensemble of n ~ 20 independent measurements. 

Unfortunately, available observational data do not allow 
us to probe AF/F at the 0.1 ppm level. First at all, only 
a handful of sources are known where both [Ci] and 13 CO 
radial velocities were measured (Schilke et al. 1995; Stark 
et al. 1996; Ikeda et al. 2002; Mookerjeaet al. 2006a,b). The 
line profiles from these observations were usually fitted with 
single Gaussians in spite of apparent asymmetries seen in 
some cases (e.g., Fig. 7 in Ikeda et al. 2002). Besides, the 
measured radial velocities were not corrected for different 
beamsizes. As a result, the scatter in AV becomes large, 
and the accuracy of the AF/F estimate deteriorates. 
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3. The a 2 /fi estimate 

In this section we consider constraints on the spatial vari- 
ations of a 2 1 \i which can be obtained from observations 
of emission lines of atomic carbon and carbon monoxide in 
submm- and mm- wave regions. The FS [C i] lines and low- J 
rotational lines of 13 CO are observed towards many galactic 
and extragalactic objects (Bayet et al. 2006; Omont 2007). 
For our purpose we selected a few molecular clouds located 
at different galactocentric distances where the radial veloc- 
ities of these species were measured with a sufficiently high 
precision (e v ~ 100 m s _1 ). 

Table [1] lists molecular clouds with both [Ci] and 13 CO 
line measurements which are available in literature. The 
data were obtained under the following conditions. 

TMC-1 — the Taurus Molecular Cloud (D ~ 140 pc). 
This dark molecular cloud was studied with the Caltcch 
10.4m submillimeter telescope on Mauna Kea, Hawaii 
(Schilke et al. 1995). The beamsize at the [Ci] (1-0) fre- 
quency was 15", while at the 13 CO (2-1) frequency it was 
about 30". Schilke et al. observed similar shapes of the [Ci] 
(1-0) and 13 CO (2-1) profiles at five positions perpendicular 
to the molecular ridge close to the cyanopolyyne peak. The 
line parameters listed in Table Q] were derived by Gaussian 
fits, although the line shapes were not exactly Gaussians. 
Therefore the errors of the line parameters are the formal 
la errors of the fitting procedure. 

LI 83 — is an isolated quiescent dark cloud at a dis- 
tance of about 100 pc (Mattila 1979; Franco 1989). The 
observations of the [Ci] and 13 CO lines at six positions 
along an east-west strip through the center of the cloud 
were obtained with the 15m James Clerk Maxwell Telescope 
(JCMT) on Mauna Kea, Hawaii (Stark et al. 1996). The 
beamsize at 492 GHz was 10" and 22" (A-band) and 15" (B- 
band) at 220 GHz. The [C i] and 13 CO data were smoothed 
to a resolution of 0.4 km s _1 and 0.2 km s _1 , respec- 
tively. These emission lines show similar asymmetric pro- 
files which can be attributed to two kincmatically different 
components closely spaced in velocity with central veloci- 
ties around 1 km s _1 and 2 km s _1 . These components are 
marginally resolved in the [C i] spectra at two positions (# 
7 and 8 in Table [l}. But since 13 CO lines were not resolved 
at these positions, we include in Tabled] the results of one 
component Gaussian fits of both 13 CO (2-1) and [Ci] (1-0) 
spectra from Stark et al. (1996). 

Ceph B — is a giant Cepheus molecular cloud at a dis- 
tance of ~ 730 pc located to the south of the Cepheus OB3 
association of early-type stars (Blaauw 1964). Cepheus B, 
the hottest 12 CO component of this complex (Sargent 1977, 
1979), is surrounded by an ionization front driven by the 
UV radiation from the brightest members of the OB3 as- 
sociation (Felli et al. 1978). The observations of the [Ci] 
(1-0) line were obtained using the KOSMA 3m submillime- 
ter telescope on Gornergrat, Switzerlaand (Mookerjea et 
al. 2006a). This dataset was complemented with 13 CO ob- 
served with the IRAM 30m telescope (Ungerechts et al. 
2000). All data were smoothed to the spatial resolution of 
1' and the velocity resolution of 0.8 km s _1 . Table Q] in- 
cludes [Ci] and 13 CO (2-1) lines arising around Vlsr of 
— 13.8 km s" 1 at the position of the hotspot in Cepheus 
B. The VLsr values of the [Ci] (1-0) and 13 CO (2-1) posi- 
tions derived from Gaussian fitting were reported in Table 2 
of Mookerjea et al. (2006a) without their errors. However, 
since the lines look symmetric (Fig. 3, Mookerjea et al. 



2006a), we assign them an error of 0.1 km s^ 1 . This is 
slightly larger than the uncertainty of ~l/10t/i of the res- 
olution element, - a typical error of the line position for 
symmetric profiles, - but does not affect significantly the 
sample mean value of AV^. 

Orion A,B — are giant molecular clouds located at 
~ 450 pc (Genzel & Stutzki 1989). The observations of the 
[Ci] (1-0) line towards 9 deg 2 area of the Orion A cloud 
and 6 deg 2 area of the Orion B cloud with a grid spacing 
of 3' were carried out with the 1.2m Mount Fuji submil- 
limeter telescope (Ikeda et al. 2002). These observations 
were complemented with the 13 CO (1-0) dataset presented 
in Table 3 in Ikeda et al.. At the frequency 492 GHz the 
spatial and velocity resolutions were, respectively, 2.2' and 
1.0 km s^ 1 , whereas at the frequency 110 GHz they were 
1.6' and 0.3 km s" 1 . The profiles of the [Ci] (1-0) and 
13 CO (1-0) lines were found to be very similar. All spectra 
were well fitted with one or two Gaussian functions, and the 
velocity centers of the [Ci] and 13 CO lines are almost the 
same: | AV | = 0.2 ±0.1 km s" 1 . The results of the Gaussian 
fitting are given in Table [T] 

Cas A — is a supernova remnant at a distance of ~ 3 
kpc (Braun et al. 1987). It was mapped in the [Ci] (1- 
0) line on the KOSMA 3m submillimeter telescope with 
the beamwidth of 55" and the velocity resolution of 0.6 
km s _1 (Mookerjea et al. 2006b). These observations have 
been compared with the 13 CO (1-0) observations (beamsize 
~ 60", spectral resolution ~ 0.1 km s -1 ) taken from Liszt 
& Lucas (1999). Both the [Ci] (1-0) and 13 CO (1-0) emis- 
sion spectra were averaged over the disk of Cassiopeia A. 
The results of Gaussian fitting of subcomponents resolved 
in the [Ci] (1-0) and 13 CO (1-0) spectra are included in 
Table [TJ Two strong emission feature observed in both [C i] 
and 13 CO (1-0) lines were identified with the Perseus arm 
at —47 km s _1 (~ 2 kpc distant) and with the local Orion 
arm at —1 km s _1 (~ 460 pc distant). 

The velocity offsets AV between the 13 CO and [C i] lines 
are given in Col. 7 of TableQ] the corresponding linewidths, 
<7„, are shown in Cols. 4 and 6. When both transitions trace 
the same material, the lighter element C should always have 
larger linewidth. If the line broadening is caused by thermal 
and turbulent motions, i.e., tr 2 = cr 2 horm + cr 2 urb , then for 
two species with masses m\ < ni2 we have 

\J m 1 /m 2 < (T2,v/<7l,v < 1 • (6) 

In practice, this inequality is fulfilled only approximately. 
Except for the pure thermal and turbulent broadening there 
are many other mechanisms which can give rise to the 
broadening of atomic and molecular lines. These are sat- 
uration broadening (lines have different optical depths), 
the presence of unresolved velocity gradients (nonthermal 
distribution is not normal), the increasing velocity disper- 
sion of the nonthermal component with increasing map size 
(the higher angular resolution is realized for the higher fre- 
quency transitions), etc. Thus, the consistency of the ap- 
parent linewidths defined by Eq. (|6j) is a necessary condition 
for two species with different masses to be co-spatially dis- 
tributed, but is not a sufficient one. 

From Table [T] it is seen that the inequality (|6]) is ful- 
filled for all selected pairs 13 CO/[Ci] within the estimated 
uncertainties of the linewidths. Thus, the whole sample of 
n = 25 AV values can be used to estimate AF/F . 

The averaging of the velocity offsets over the dataset 
gives the unweighted mean AV^ W = (Vls_r( 13 CO) — 
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Table 1. Parameters derived from Gaussian fits to the 13 CO J = 2 — 1, J = 1 — and [Ci] J = 1 — emission line 
profiles observed towards Galactic molecular clouds. V LSR is the line center, a v is the line width (FWHM). The numbers 
in parentheses are the standard deviations in units of the last significant digit (see text for more details). Col. 7 lists 
velocity offsets = V°° R — V°/ R , and their estimated errors. 



Source 


No. 


yCO 


<?v > 


v LSRl 


„ c 1 


AV, 


Ref. 






km s _1 


km s _1 


km s _1 


km s _1 


km s 1 




<i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



TMC-1 



L183 



Cas A 



1 


6.K 


T 


.a 


2.01 


1) 


2 


6.1| 


T 


.a 


1.61 


"1 


3 


6.11 


1* 


.a 


1.71 


"lj 


4 


6.11 


T 
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Vlsr(Ci)) = 0.046 ± 0.083 sta t km s" 1 . With weights 
inversionally proportional to the variances, one derives 
AV W = 0.029 ± 0.077 km s _1 . The median of the sam- 
ple is AVmod = 0.0 km s _1 , and the robust M-estimate 
(L10) is AV M = 0.022 ± 0.082 km s" 1 . The statistical 
error for the mean velocity offset measurement is larger 
than that expected from the published values of the sta- 
tistical errors from the one component Gaussian fits: the 
mean error of the individual AV is 0.13 km s -1 , and the 
expected error of the mean AV is ^0.026 km s _1 . A pos- 
sible reason for such a high Doppler noise has been dis- 
cussed in Sect. 2. The systematic error in this case is dom- 
inated by the uncertainty of the rest frequency of the [C 1] 
(1-0) transition, e v = 33.5 m s _1 . Thus, taking the M- 
cstimate as the best measure of the velocity offset, we have 
AV M = 0.022 ± 0.082 sta t ± 0.034 sys km s" 1 , and the la 
upper limit on |AV| < 0.11 km s _1 . 

This estimate restricts the spatial variability of F at 
the level of \AF/F\ < 0.37 ppm. Recently we obtained a 
constraint on the spatial change of the electron-to-proton 
mass ratio A/i//x < 0.03 ppm based on measurements in 
cold molecular cores in the Milky Way (L10). Combining 
these two upper limits, the fine-structure constant can be 
bound as |Aa/a| < 0.2 ppm. 

4. Conclusion 

The level of 0.2 ppm represents a model-dependent up- 
per limit on the spatial variations of a. Under model- 
dependence we assume here that both AF/F and Aji/ji 



do not change significantly from cloud to cloud, since as- 
trophysical measurements of these parameters are made 
in low density regions of the interstellar medium with 

Pcloud ^ Pterrestrial- 

For comparison, the upper limit on the temporal co- 
variation obtained from high-redshift quasar absorbers is 
|Aa/a| < 2 ppm (Sect. 1). If dependence of constants on 
the ambient matter density dominates over temporal (cos- 
mological) , as suggested in chameleon-like scalar field mod- 
els, then one may expect that |Aa/a| < 0.2 ppm at high 
rcdshifts as well, since quasar absorbers have gas densi- 
ties similar to those in the interstellar clouds. Taking into 
account that the predicted changes in a and \x are not inde- 
pendent and that ^.-variations may exceed variations in a 
(e.g., Calmet & Fritzsch 2002; Langacker et al. 2002; Dine 
et al. 2003; Flambaum et al. 2004), even a lower bound of 
\Aa/a\ < 0.03 ppm is conceivable within the framework of 
the chameleon models. 

We note that if a theoretical prediction \Aa/a\ 
|A/x//x| is valid, then AF/F w —A/j,/(j,, and, hence, the F- 
cstimate with a further order of magnitude improvement in 
sensitivity will provide an independent test of the tentative 
change of /x. 

The factors limiting accuracy of the current estimate of 
AF/F at z = are a relatively low spectral resolution of 
the available observations in submm- and mm- wave bands, 
a rather large uncertainty of the rest frequencies of the [C 1] 
FS lines, and a small number of objects observed in both 
[Ci] and 13 CO transitions. 
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Modern telescopes like the recently launched Herschel 
Space Observatory can provide for Galactic objects the 
spectral resolution as high as 30 m s _1 (e.g., the 
Heterodyne Instrument for the Far Infrared, HIFI, has 
resolving power R = 10 7 ). This means that the posi- 
tions of the [Ci] FS lines can be measured with the 
uncertainty of ^3 m s _1 . In the near future, high 
precision measurements will be also available with the 
Atacama Large Millimeter/submillimeter Array (ALMA), 
the Stratospheric Observatory For Infrared Astronomy 
(SOFIA), the Cornell Caltech Atacama Telescope (CCAT) 
and others. Thus, any further advances in exploring AF/F 
depend crucially on new laboratory measurements of the 
[C i] FS frequencies. If these frequencies will be known with 
uncertainties of a few m s _1 , then the parameter AF/F 
can be probed at the level of 10~ 8 which would be com- 
parable with the non-zero signal in the spatial variation of 
the electron-to-proton mass ratio /i. 
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